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Abstract 

The singularity structure and the corresponding asymptotic behavior of a 3-brane 
coupled to a scalar field in a five dimensional bulk is analyzed in full generality, 
using the method of asymptotic splittings. It is shown that the collapse singularity 
at a finite distance from the brane can be avoided only at the expense of making 
the brane world-volume positively or negatively curved. 
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1 Introduction 



Some time ago, an interesting idea to address the cosmological constant problem was 
proposed, based on the so-called self-tuning mechanism [TJ |2]. The simplest model con- 
sists of a 3-brane embedded in a five- dimensional bulk, in the presence of a scalar field. 
The later is coupled to the brane in a particular way, motivated by string theory, that 
allows flat brane world-volume solution independently of the brane tension value. It was 
however realized that a singularity appears in the bulk, at some finite distance from the 
brane, which can also be thought as a reservoir through which the vacuum energy decays. 

An obvious question is then whether the development of such a singularity is a generic 
feature of these models, or under what conditions may be avoided. Here and in two 
subsequent papers, we investigate this question in a generalized class of models. Since 
in this case a general solution cannot be found analytically, we use a powerful tool 
developed a few years ago, called method of asymptotic splittings, that allows to compute 
all possible asymptotic behaviors of the equations of motion around the assumed location 
of a singularity [3J. Our main result is twofold: 

• The existence of a singularity at a finite distance is unavoidable in all solutions 
with a flat brane. This confirms and extends the results of earlier works that made 
similar investigations in different models, using other methods [U |5]. 

• The singularity can be avoided (e.g. moved at infinite distance) in several cases 
where the brane becomes curved, either positively or negatively. Thus, requiring 
absence of singularity brings back the cosmological constant problem, since the 
brane curvature depends on its tension that receives quartically divergent quantum 
corrections. 

As mentioned already, our results are established in an extended version of the simplest 
model, where the bulk scalar field has a general coupling to the brane, motivated for 
instance from a loop corrected string effective action, allowing for curved world-volume. 
Besides the scalar field (or in the absence of it), we consider a generic bulk matter content 
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parametrized by a fluid with an arbitrary equation of state. For pedagogical purposes, 
our analysis is separated in three parts contained in three different publications. A 
preliminary version of our results was published in [6]. 

In this first paper, we give a detailed picture of the dynamical evolution of an extended 
version of the simplest model to include curved branes. We show that the emergence of 
the finite-distance singularity is the only possible asymptotic behavior for a flat brane, 
whereas for a curved brane the singularity is shifted at an infinite distance. We provide a 
detailed study of the asymptotics of this model using the method of asymptotic splittings 
expounded in [3]. 

The structure of this paper is as follows: In Section 2, we derive the form of the 
dynamical system on which our subsequent asymptotic analysis is applied. In Sections 
3 and 4, we give the asymptotics of the models consisting of flat and curved brane 
respectively. In Section 5 we conclude and also comment on possible future work in 
various directions, considering for instance other forms of matter in the bulk. In the 
Appendix, we briefly outline the basic steps of the method of asymptotic splittings. 

2 Dynamics of scalar field-brane configuration 

In this Section we set up the basic equations for our braneworld. We study a braneworld 
model consisting of a three-brane embedded in a five- dimensional bulk space with a scalar 
field minimally coupled to the bulk. The total action Stotai splits in two parts, namely, 
the bulk action S^uik and the brane action Sfc rane , 



total Sbulk 'Sferanej 



(2.1) 



with 




(2.2) 



(2.3) 
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where Y denotes the fifth bulk dimension, Y* is the assumed initial position of the brane, 
A is a parameter defining the type of scalar field 0, k\ = M^ 3 , M 5 being the five- 
dimensional Planck mass, and /(0) denotes the tension of the brane as a function of the 
scalar field. 

Varying the total action (12. ip with respect to g , we find the five-dimensional Ein- 
stein field equations in the form [6], 

R AB - \u B R = X4 (v.0V B - \ 9AB (V4,f) + -^J^SM SV 

while the scalar field equation is obtained by variation of the action (12. ip with respect 
to [6] and it is: 

AD 5 = j== — 6{Y), (2.5) 

Vdet# 5 6(p 

where A,B = 1,2,3,4,5 and a, (3 — 1,2,3,4 while S(Y) = 1 at Y = Y* and vanishing 
everywhere else, and 

v A (^^g AB v B <p). (2.6) 



(2.4) 



A/det^ 

In the following we assume a bulk metric of the form 

g 5 = a 2 (Y) 9i + dY 2 , (2.7) 
where is the four dimensional flat, de Sitter or anti de Sitter metric, i.e., 

g 4 = -dt 2 + f 2 g 3 , (2.8) 

where 

33 = dr 2 + h\g 2 (2.9) 

and 

g 2 = d9 2 + sin 2 6d<p 2 . (2.10) 

Here fk = l,cosh(Ht)/H,cos(Ht)/H (if -1 is the de Sitter curvature radius) and hk = 
r, sin r, sinh r, respectively. 
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The field equations Il2.4l) - fl2.5p then take the form 

4 = (2.U) 
a z 12 a 2 

a - = (2.12) 
a 4 

0" + 4-0' = 0, (2.13) 
a 

where a prime denotes differentiation with respect to Y and k = 0, ±1. The variables 
to be determined are a, a' and <$' . These three equations are not independent since Eq. 
(I2.12p was derived after substitution of Eq. (I2.1ip in the field equation G aa = nlT aa , 
a = 1,2,3,4, 

a" a' 2 kH 2 ~<j)' 2 , N 

- + — - — = -Ak l^-. 2.14 
a a 2 a 1 6 

In our analysis below we use the independent equations ( 12. 12ft and 1 12. 13p to determine 
the unknown variables a, a' and <f)', while Eq. (12. lift will then play the role of a constraint 
equation for our system. 

Assuming a Y — > —Y symmetry and solving the Eqs. (12.41) (the -aa- component, 
a = 1, 2, 3, 4) and ( 12. 5 j) on the brane we get 

o'(y.) = -^}(<t,(Y,))a(Y.), (2.15) 

W) = M. (2.16) 

The particular coupling used in [1] allows only for flat solutions to exist. This easily 
follows by using equations ( 12. 15ft and 1 12. 16ft and solving the FRW equation ( 12. lip on the 
brane for kH 2 : 

Clearly, k is identically zero if and only if: 

7W = 2 V IT 5 - 

or equivalently, if and only if /(</>) oc e 21 /^^ 5 ^ (the authors of p] have set A = 3 and 
hence the appropriate choice for the brane tension in that case is /(</>) oc e 2ft5( ^). In our 
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more general problem, the coupling function cannot be fixed this way. By working with 
other couplings we can allow for non-flat, maximally symmetric solutions to exist and 
avoid having the singularity at a finite distance away from the position of the brane. 

For the rest of this paper our purpose is to find all possible asymptotic behaviours 
around the assumed position of a singularity, denoted by Y s , emerging from general or 
particular solutions of the system (I2.lip -f l2.13p . The most useful tool for this analysis 
is the method of asymptotic splittings [3] (see the Appendix for a brief introduction) in 
which we start by setting 

x = a, y = a', z = 0', (2-17) 

The field equations (12.121) and (I2.13P become the following system of ordinary differential 
equations: 

x = y (2.18) 
y = -\Az 2 x (2.19) 



-Ay-, (2.20) 

x 



where A = nf/A. Hence, we have a dynamical system determined by the non-polynomial 
vector field 

f = (y,-\Az 2 x,-4y^y . (2.21) 

Equation (12. lip does not include any terms containing derivatives with respect to T; it 
is a constraint equation which in terms of the new variables takes the form 



y 2 AX 2 kH 2 



2 , z + -ZT- ( 2 - 22 ) 



Equations (12. 18j) - f f2T20|) and (12.221) constitute the basic dynamical system of our study. 
There are two major cases to be treated, the first is when we choose k = in (12.221) 
and corresponds to a brane being flat, while in the second case k ^ 0, giving constant 
curvature to the brane. We shall treat these two cases independently in what follows. 
One important result of our analysis of this system will be that the inclusion of nonzero 
curvature for the brane moves the singularity an infinite distance away from the brane. 
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3 Flat brane: Finite-distance singularity 



In this Section we take k = in the basic constraint equation 



(3.1) 



i' 3 

We shall show that the only possible asymptotic behaviour of the solutions of this system 
(flat brane) is that a — > 0, a' — > oo and <p' — > oo, as Y — > Y s . 
We start our asymptotic analysis by inserting the forms 



(x,y,z) = {aV,pT>,6'r), 



in the system (l2TTgj) - (l230j) . where 

(p,q,r)eQ 3 and (a, (3, 5) E C 3 \ {0}. 



(3.2) 



(3.3) 



We find that the only possible dominant balance in the neighborhood of the singularity 
(that is pairs of the form B = {a, p}, where a = (a, (3,5), p = {p, q, r), determining the 
dominant asymptotics as we approach the singularity) is the following balance 



B 1 = {(a, a/A, V3/(aVAX)), (1/4, -3/4, -1)}. 



(3.4) 



(A second balance £>2 becomes only possible when we allow for non-zero curvature, k ^ 0, 
and will be analysed in the next Section. There are no other acceptable balances, hence 
all the possible asymptotic behaviours for a flat and curved brane can be described 
uniquely by the balances B\ and B 2 respectively.) 

Let us now focus on building a series expansion in the neighborhood of the singular- 
ity to justify the asymptotics found above. We start by calculating the Kowalevskaya 
exponents, eigenvalues of the matrix /C = Df (a) — diag(p), where Df (a) is the Jacobian 
matrix of f , which in our case reads: 

/ n i n \ 



Df(x,y,z) 





-A\z 2 
Ayz 







-2A\xz 

x x 



(3.5) 
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to be evaluated on a. For the B\ balance we have that a 
p = (1/4,-3/4,-1), thus 

/ 1 



K 



4 
3 

~16 

V3 



1 

3 
4 







a 



a, a/4, a/3/(4a/AX)), and 
\ 

(3.6) 



V 4av^4A 
The /C-exponents are then given by 



aVAX 



spec(/C) = {-1,0,3/2}. 



(3.7) 



These exponents correspond to the indices of the series coefficients where arbitrary con- 
stants first appear. The —1 exponent signals the arbitrary position of the singularity, Y s . 
We see that the first balance B\ has two non-negative rational eigenvalues which means 
that it describes the asymptotics of a general solution in the form of a series expansion, 
i.e., a series form of the solution having the full number of arbitrary constants (which for 
our system equals to two) . In order to construct an asymptotic expansion of this solution 
valid in the neighborhood of the singularity, we substitute in the system fl2.18p - fl2.20p 
the series expansions 

x = T p (a + S^ lCj T^), 

where x = (x, y, z), Cj = (c-,i, c.,-2, c^), and s is the least common multiple of the denom- 
inators of the positive eigenvalues; here s = 2, and the corresponding series expansions 
are given by the following forms: 



x 



V°o W2+1/4 ,._yoo T j/2-3/4 



y°° rJ J '/ 2_1 



Therefore we arrive at the following asymptotic solution around the singularity: 

x = aT 1 ^ + ! C32 T 7 /4 + . . . 
y = ^T~ 3/4 + c 32 T 3/4 + • - • 



-1 



4^3 
la\[A\ 



C32 



T l/2 + 



(3.8) 

(3.9) 
(3.10) 

(3.11) 
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The last step is to check whether for each j satisfying j/2 = p with p a positive 
eigenvalue, the corresponding eigenvector v of the K. matrix is such that the compatibility 
conditions hold, namely, 

(3.12) 



v 1 ■ Pj = 0, 



where Pj are polynomials in Cj, . . . , Cj_i given by 

/Cc, - (j/s) Cj = P y (3.13) 

Here the relation j/2 = 3/2 is valid only for j = 3 and the associated eigenvector is 

ay/AX 7a\/AX 



v 



The compatibility condition, 



y/3 ' AV3 



v{K- (3/2)X 3 )c 3 = 0, 



(3.14) 



(3.15) 



therefore indeed holds since 

/ 7 



(K - (3/2)X 3 )c 3 = c 32 



\ 



_ 3 
~16 



3 
4 



a 



V3AX 



3 
2 



/ 



4 
7 



4v/3 



o 



V / 

(3.16) 

This shows that a representation of the solution asymptotically by a Puiseux series as 
given in Eqs. (13.9j) -( l37Tl~]) is valid. Hence we conclude that near the singularity at finite 
distance Y s from the brane the asymptotic forms of the variables are: 



a 0, a' — > oo, 0' — > oo. 



(3.17) 



This is exactly the asymptotic behaviour of the solution found previously by Arkani- 
Hammed et al in pQ. Our analysis shows that this is the only possible asymptotic 
behaviour for a flat brane since there exist no other dominant balances in this case. 
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4 Curved brane: Infinite-distance singularity 



In this Section we show that the collapse singularity that necessarily arises in the case 
of a flat brane is avoided (or shifted at an infinite distance away from the brane) when 
we consider a curved brane instead. 

The new asymptotics follow from the study of a second balance that results from the 
substitution of ( 13. 2 p in (I2.18p - ()2.20p . We calculate this new balance to be, 

B 2 = {(a,a,0),(l,0,-1)}. (4.1) 

It corresponds to a particular solution for a curved brane since it satisfies Eq. ( 12. 22ft for 
k and a 2 = kH 2 (here we have to sacrifice one arbitrary constant by setting it equal 
to kH 2 ), k — ±1. The JC- matrix of B2 is 





/-I 


1 


\ 


K = Df((a, a, 0)) - diag(l, 0, -1) = 













V 





"3/ 



with eigenvalues 

spec(/C) = {-1,0,-3}. (4.3) 

Thus for the balance £>2 we find two distinct negative integer /C-exponents and an infinite 
expansion in negative powers of a particular solution (recall that we had to sacrifice 
one arbitrary constant) around the presumed singularity at Y s , with the negative K,- 
exponents signaling the positions where the arbitrary constants first appear [7]. We 
therefore expand the variables in series with descending powers of T in order to meet 
the two arbitrary constants occurring for j = —1 and j = —3, i.e., 

x = S7~ Cil T^ 1 , y = E i -~c, 2 T', z = EJ^sT'" 1 . (4.4) 

Substituting these series expansions back in the system (I2.18p - (l2.20p and after some 
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manipulation, we find the following asymptotic behaviour, 



x = aT + c_i i + 
y = a H 

z 



C- 33 T- 4 + 



(4.5) 
(4.6) 
(4.7) 



Let us check the compatibility conditions for j — — 1 and j 

( 1 W c_u \ 

(/C+Z 3 )c_ 1 = 



and 



1 
\0 -2 J 

( 2 1 \ 



(/C + 3X 3 )c_ 





V / 
/ o \ 



-3. We find that 
, 





V C 33 / 



( \ 



w 



(4.9) 



3 
\ / 

so that the compatibility conditions are indeed satisfied. The expansions given by Eqs. 
( 14.5p -( l4~7l) are therefore valid, and we can say that as T — > 0, or equivalently as S = 
1/T — > oo, we have that 



a — > oo, a — > oo, 



— > oo. 



(4.10) 



Therefore for a curved brane we find that there can be no finite-distance singularities. 
The only possible asymptotic behaviour is the one given in f)4.10p which is only valid at 
an infinite distance from the brane. 



5 Conclusions 

In this paper we studied a braneworld consisting of a three-brane embedded in a five- 
dimensional bulk space filled with a scalar field with a special emphasis in the possible 
formation of finite-distance singularities away from the brane into the bulk. We have 
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shown that the dynamical behaviour of this model strongly depends on the spatial ge- 
ometry of the brane, in particular whether it is flat or not. For a flat brane the model 
experiences a finite-distance singularity toward which all the vacuum energy decays (since 
<p' —7- oo, as Y — > Y s ), whereas for a curved brane the model avoids the singularity which 
is now located at an infinite distance. 

It is interesting that a third balance which initially arises from the substitution of 
(E2]) in <KW\ - <K2H\f . namely, the form 

£ 3 = {(«, 0,0), (0,-1,-1)}, 

is not acceptable for the model we consider in this paper since it does not give the nec- 
essary — 1 /C-exponent. In future work (8j and [9], we will see that this balance although 
impossible in the case treated here, does become possible (although in a somewhat 'mild' 
form) when we replace the scalar field studied here with other matter components such 
as a perfect fluid or a combination of a perfect fluid and a scalar field. We therefore 
conclude that for the case of interest in this paper the collapse singularity found is the 
only type of singularity that can develop at a finite distance from a flat brane. 
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A Appendix: The method of asymptotic splittings 

We refer briefly here to the basic steps of the method of asymptotic splittings. A detailed 
analysis can be found in Ref. [3]. 
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Consider a system of n first order ordinary differential equations 

x = f(x), (a.i; 



where x = (xi, . . . ,x n ) E MJ 1 , f(x) = (/i(x), . . . , / n (x)) and ' = -4?, Y being the inde- 
pendent variable. In this paper, we refrain from calling Y a time variable and giving it 
the interpretation of time. Since we are interested in singularities located at a distance 
from the brane and into the bulk, it seems more appropriate to talk about finite-distance 
singularities and give to the Y variable a spatial interpretation. The general solution of 
the above system contains n arbitrary constants and describes all possible behaviours of 
the system starting from arbitrary initial data. Any particular solution of ( lA.lj) . on the 
other hand, contains less than n arbitrary constants and describes a possible behaviour 
of the system emerging from a proper subset of initial data space. 

We say that a solution of the dynamical system (1A.ll) exhibits a finite-distance sin- 
gularity if there exists a Y s e R and a xo G W 1 such that 

lim ||x(F;xo)|| ^00, (A.2) 

where . is any LP norm. The purpose of singularity analysis (cf. [3], [10]) is to build 
series expansions of solutions around the presumed position of a singularity at Y s in order 
to study the different asymptotic behaviours of the solutions of the system ( lA.lj) as one 
approaches this singularity. In particular, we look for series expansions of solutions that 
take the form of a Puiseux series (any log terms absent), namely, a series of the form 

x = T p (a+S~ 1 c i T'/ s ), (A.3) 

where T = Y - Y s , p e Q n , s G N. 

The method of asymptotic splittings for any system of the form (lA.lj) is realized by 
taking the following steps: 

• First, we find all the possible weight-homogeneous decompositions of the vector field 
f by splitting it into components f^: 

f = f(°) + f « + ... + f( fc ), (A.4) 
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with each of these components being weight homogeneous, that is to say 

f (j) (aT p ) = T p+1( - qU) - l) { u) (8L) j = 0,...,k, (A.5) 

where ael" and q*' are the positive non-dominant exponents that are defined by (I A. 71) 
below. 

• We substitute the forms x = aT p in the system x' = f^(x) in order to find all 
possible dominant balances, i.e., finite sets of the form {a, p}. The order of each balance 
is defined as the number of the non-zero components of a. 

• For each of these balances we check the validity of the following dominance condi- 
tion ' 

fW(aTP) 

lim J ~* \ = 0, A.6) 

and define the non-dominant exponents j = 1, . . . , k by the requirement that 

T qU) . (A. 7) 



E fe =if 0-)( T P) 



Tp- 1 

The balances that cannot satisfy the condition (1A.6I) are then discarded. 

• We compute the Kovalevskaya matrix /C defined by 

K = L>f (0) (a) -diagp, (A.8) 

where Df^(a) is the Jacobian matrix of evaluated at a. 

• We calculate the spectrum of the /C-matrix, spec(fC), that is the set of its n eigen- 
values also called the IC-exponents. The arbitrary constants of any particular or general 
solution first appear in those terms in the series (1A.3I) whose coefficients have indices 
k = ps, where p is a non-negative /C-exponent and s is the least common multiple of the 
denominators of the set consisting of the non-dominant exponents q^ and of the positive 
/C-exponents (cf. [3], [I0])- The number of non-negative /C-exponents equals therefore 
the number of arbitrary constants that appear in the series expansions of (IA.3I) . There 
is always the —1 exponent that corresponds to the position of the singularity, Y s . (A 
dominant balance corresponds thus to a general solution if it possesses n — 1 non- negative 
/C-exponents (the nth arbitrary constant is the position of the singularity, Y s )). 
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• We substitute the Puiseux series: 

x i = E^ c ii T»^>i = l,...,n ) (A.9) 

in the system (1A.1I) . 

• We find the coefficients Cj by solving the recursion relations 

JCcj - j -cj = P j (c 1 ,..., Cj_i) (A. 10) 

where Pj are polynomials that are read off from the original system. 

• We verify that for every j = ps, with p a positive /C-exponent, the following 
compatibility conditions hold: 

v 1 -P j = 0, (A.ll) 

where v is an eigenvector associated with the positive /C-exponent p. 

• We repeat the procedure for each possible decomposition. 

We note that if the compatibility condition above (Eq. (1A.11I) ) is violated at some 
eigenvalue in the spec(/C), then the original Puiseux series representation of the solution 
cannot be admitted and instead we have to use a ip-series for each one of the eigenvalues 
with this property. This is a series that includes log terms of the form 

x = T p (a + E^. =lCii T i/s (T"logT) J >) , (A.12) 

where p is the /C-exponent for which the compatibility condition is violated. The rest of 
the procedure in this case is the same as before. 
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